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Abstract 

We calculate cross sections of high energy electron inclusive scattering off 
nuclear matter in a new and consistent formulation based on the Green's function 
method with the Glauber approximation, which is an extension of our previous 
work on the nuclear transparency in (e, e'p) reaction. The comparison with other 
approaches is discussed. In this framework, we study the finite-range effect of the 
nucleon-nucleon interaction in the final-state interactions, and the effect of the 
nuclear short-range correlation. We propose a zero-range approximation, which 
works well when correlation and finite interaction-range effects are included. It 
greatly reduces the numerical work, while maintaining a reasonable accuracy. 
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1 Introduction 



Experimental data on high-energy electron-nucleus inclusive scattering have been ac- 
cumulated for a decade Q. The experiments were mostly done at SLAC, and will be 
done at the Jefferson Lab. Also a new facility called MUSES (Multi-Use Experimental 
Storage ring) is planned at RIKEN ||. Precise data will become available with these 
facilities and provide us with a way of studying where and how the quark-gluon degrees 
of freedom may come into sight. 

The process that we study is the (e, e') inclusive reaction (Fig. 1) in the nucleon 
quasi-elastic region, i.e., 

e + e' + N + (A-iy. (1.1) 

The incident electron has a four-momentum, = (e, k e ) , and the outgoing electron, k'^ 
= (e', k' e ). The momentum transfer is — — k' = (u, q). We are interested in the 
quasi-elastic region, i.e., |q| ~ 2 [GeV/c], and to < 2 [GeV]. So the elementary process 
is considered to be mostly the electron-nucleon elastic scattering, e + N —> e + N, and 
the final-state interaction (FSI) of the struck nucleon is what we wish to study here. 

In this work we calculate the effect of the FSI, taking into account only the nucleonic 
degrees of freedom. We use the Glauber approximation for the FSI, including the 
nuclear correlation. We believe that our treatment is a more systematic approach than 
the previous works. 

The theoretical treatment of this reaction is similar to that we used for the nuclear 
transparency in the (e,e'p) reaction at large-momentum transfers 0, [|]- This process 
has recently attracted many researchers in connection with a speculated phenomenon 
of the color transparency [[J. The present formulation allows us to describe the (e, e') 
and (e, e'p) reactions in a unified way. 

The nuclear transparency for electron scattering can be defined as the ratio of 
the (e, e'p) response to the inclusive (e, e') response in the quasi-elastic region, and is 
conveniently formulated by the Green's function method We will give here a way 
of calculating the response functions based on the Glauber approximation for the 
Green's function, without taking account of the proton internal structure. 

The contents of this paper are as follows: In sec. |2] we explain our formulation. We 
explain our standpoint about the elementary process, i.e., the electron-nucleon cross 
section, in subsec. pTT[ We derive a closed form of the inclusive cross section based on 



the Glauber approximation in subsec. 2.2. The definition of the FSI-function is also 



given there. The expression of the response function for nuclear matter is derived in 
subsec. 2JB, and the definition of our "convolution" function is found. We derive the 
approximate expressions of the FSI-function when we apply zero-range approximations 
for the nucleon-nucleon potential in subsec. [2.4| . In sec. |3] we compare our treatment of 
the FSI with others. In subsec. [3.1| we discuss the difference and the relation between 
our formulation and the theory of Gersch, Rodriguez, and Smith M, M, [f[0[|. A similar 



comparison with the optical potential formalism [11], [12 is given in subsec. 3^2. In 



sec. 0] we show our numerical results and discuss their indications. We show the FSI- 



function in subsec. 4.1, and the inclusive cross sections in subsec. 4.2. We summarize 



our results in sec. |5|. In the Appendix we show a method for deriving a nucleon-nucleon 
potential from the experimental scattering amplitude in the eikonal approximation. 
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2 Formulation 



In this section we explain our formulation of the electron-nucleus inclusive scattering 
based on the Green's function method 0. We employ the Glauber approximation for 
the FSI of the struck nucleon. 



2.1 Choice of the eN Cross Section 

Let us first explain our standpoint about the elementary process of the (e, e') inclusive 
reaction. 

We start with the (e, e') inclusive differential cross section on a target nucleus, A, 
which is assumed to be expressed as [ p~3|] 

da eA ,da eN 

(-jr)ei,on%q), (2.1) 



dVtduo dVt 

where the cross section on the r.h.s. is that of on-shell electron-nucleon elastic scat- 
tering for the same incident energy and scattering angle, and (...) implies averaging 
over the spin, and the isospin, of the target nucleons. The difference between the 
longitudinal and the transverse responses is also neglected. 

This factorized form, eq. ( |2.1|) , is an approximation, on which some comments are 
in order: The nucleon struck by the electron is not in free space, but in the nuclear 
medium both in the initial and in the final states. Such a nucleon is often said to be 
off-the-mass shell, and its electromagnetic form factors to be used in the calculation of 
the cross section are generally different from those of the on-shell nucleon. However, 
we have no reliable way of estimating the differences. In fact, we need a model for 
the internal structure of the nucleon as well as a model of its interaction with the 
nuclear medium to study the off-shell effects. We note here that the off-shell effects are 



essentially dynamical and cannot be obtained by kinematical considerations [[14]], (15 
In the following we discuss how off-shell the nucleon will be in the case of inclusive 
(e, e') on nuclear matter. 

Let us assume that the single nucleon spectrum in the nuclear matter is given by 

(E-V) 2 = (m N + S) 2 + p 2 , (2.2) 

where E, p and m-^ are the energy, the spatial momentum and the mass of the nucleon, 
respectively, with V (S) being the single particle vector (scalar) potential. The degree 
of the nucleon being off-the-mass shell is measured by 

5m 2 = E 2 - p 2 - ml = 2EV - V 2 + 2m N S + S 2 . (2.3) 



The Walecka model gives V = 300 [MeV], S = -350 [MeV] (Case I). We will also 
consider the two extreme cases, i.e., V = [MeV], S = —50 [MeV] (Case II), and 
V = -50 [MeV], S = [MeV] (Case III), keeping the sum, S + V, to be the same. 
Figure 2 shows Sm 2 /m 2 ^ as a function of E for these three cases. In all the cases, 
the nucleon in the initial state, where E ~ m-^, is not very off-shell, because 5m 2 /m 2 ^ 
~ —0.1. On the other hand, the nucleon in the final state, where E ~ 2 — 3m^, is 
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far off-shell, i.e., Sm 2 /m^ ~ 0(1), in the case (I), while it is not very off-shell in the 
cases (II) and (III). This exercise tells us that one should be careful in choosing the 
off-shell kinematics for the inclusive (e, e') response. What may be important in the 
high-energy inclusive reaction is the off-shellness of the nucleon in the final state, not 



in the initial state. This has been stressed previously by, e.g., Uchiyama, et al. [17 



and Ciofi and Simula ITS 



The off-shell cross section proposed so far is that of the half-off-shell [Q, flSfl . In 
the prescription the bound nucleon in the initial state is treated as off-the-mass shell, 
while the nucleon in the final state is treated on-the-mass shell. This is reasonable for 
the (e, e'p) semi-inclusive reaction, but not for the (e, e') inclusive. Furthermore, for the 
(e, e') inclusive process the nuclear matter limit (A — > oo) can be taken, while for the 
(e, e'p) semi-inclusive process that limit cannot be defined, because the nucleon in the 
final state cannot be free due to the infinitely extended nuclear matter. Therefore, the 
off-shell cross section, which is presently available, is not suitable for the elementary 
process of the (e, e') inclusive reaction. 

We thus choose to avoid using such off-shell cross sections, and use the on-shell cross 
section. This is the basis for our choice of the on-shell kinematics for the initial nucleon. 
The construction of the models describing the internal structure of the nucleon and its 
interaction with the nuclear medium is beyond the scope of this work, and we use the 
simple factorized form with the on-shell form factors. 

Therefore the agreement of the numerical results with the experimental data is 
not necessarily the aim of this work. In the following we concentrate on a consistent 
treatment of the FSI of the struck nucleon in the Green's function method with the 
Glauber approximation. 



2.2 Green's Function Method 

We now proceed to the derivation of the closed form of the response function, S(uj, q), 
in eq. Q2.1| ) in the Green's function method with the Glauber approximation. The 
response function is defined by 

5(w,q) = Z\(X\d(q)\A)\* 8(u> - E x ) 
x 

= --Im%q), (2.4) 

7T 

where 

R(u, q) = (A\6\q) i 0(q) \A), ( v > 0). (2.5) 

oj — h A + ir] 

The initial state, A, is the ground state of the target nucleus, and the sum over the 
final states, X, is taken, because we are considering the inclusive process. Ha is the 
full hamiltonian of the target nucleus. The energy of the final state, Ex, in eq. ( |2.4| ) 
is measured from that of the ground state of the target. 

0(q) in eq. fl2.4|) is a hard interaction operator defined by 

6(q) = exp{zq-ri}, (2.6) 
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where ri is the coordinate of the struck nucleoli denoted as 1. We neglect the interfer- 
ence effects among different nucleons with which the photon couples. 
We decompose the full hamiltonian, H A , in eq. (|2.4|) as 



H A = Ki + Ha-i + Via-i- (2.7) 

Here K\ is the kinetic-energy operator for the struck nucleon, 1, and H A _i is the 
hamiltonian for the residual nucleus. V\,a-\ is the interaction of the struck nucleon 
with other nucleons in the target nucleus, which is assumed to be expressed as a sum 
of the two-body interactions, VnnM, 

A 

^ M -iOi; r 2 , • • • , r A ) = £ ^ NN (n - r,). (2.8) 

i=2 

For numerical calculations we construct the nucleon- nucleon potential, Vnn(*"), from 
the phenomenological scattering amplitude describing the observed cross sections. The 
potential is complex in the energy region where we are interested in. We will leave the 
details in the Appendix. 

Introducing the Green's function for the nucleon, 1, interacting with the residual 
nucleus in the fixed-scatterer approximation which is given by 

G(u;r 1 ,r' 1 ;r 2 ,...,r A ) = (n| T \ -|r'i), (2.9) 

UJ-Kx- V 1)A -1 + IT) 

the response function is given as 

S(uj,q) = — — Im / dr 1 dr' 1 dr 2 ■ ■ ■ dr A exp{zq ■ (r\ — i^)} (2.10) 

7T J 

x *^(ri, r 2 , . . . , r A ) ^(r'i, r 2 , . . . , r A ) G(uj; n, r'x; r 2 , . . . , r A ), 

where \l/(r 1 ,r 2 , . . . , r^) is the A-body nuclear wave function. The fixed-scatterer ap- 
proximation implies that the excitation energy, E n , of the nucleus is neglected, and 
that the nucleon-nucleon interaction is assumed to be local [ ]1| . The struck nucleon is 
considered to become energetic enough to justify the fixed-scatterer approximation for 
the FSI due to Vnn(i")- 

The Green's function is now a one-body Green's function with a multi-centered 
potential, and we here introduce the eikonal expression for it 



G(u;r 1 ,r' 1 ;r 2 ,---,r A ) 

A 

_I £( bl - b'O B{z x - z[) exp{ip(z 1 - z[) - - r dz'l V NN (r"i - r,)} 
+ (backward piece), (2-H) 
where p is related to the energy loss, lu, by 



u + m N ~ E p = Jp 2 + mjj, (2.12) 
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with mu being the nucleoli mass, and v = p/E p . We are free to choose the direction of 
the z-axis at this stage, but will take it to be in the direction of q, since the exponential 
in the integrand of eq. fl2.10|) is strongly oscillating in this direction for large- |q|. Here 
we have also neglected E n , because it is small compared to u, which is roughly 2 [GeV]. 
The vectors, r 1; r' 1; and r"x, are thus decomposed as, 

r 1 = (bi^i), v' 1 = (b' 1 ,z' 1 ), r" 1 = (bi,^)- (2-13) 

The backward piece in eq. ( |2.11| ) becomes important only when we discuss the sum 
rule, and will be omitted in the following. Then the imaginary part of the Green's 
function is obtained as 

Im G(u; ri, r' x ; r 2 , • • • , r A ) = — {G(u; n, r\; r 2 , • • • , r A ) - G*(u; r\, n; r 2 , • • • , r A )} 
= — 5(bi - b'i) exp{zp(zi - z[)} 



% r z i 

:[9(zi - z[) exp{ — / dz" 22 Vnn(i""i - r^)} 

V Jz 'l j=2 

+ 6(z[ - Zl ) exp{- ^ dz'{ E VUr"i - r,-)}]. (2.14) 



Substituting this expression into the r.h.s. of eq. ( |2.10| ), we obtain the response func- 
tion, S(u,q), in the Glauber approximation. 

Although the response function is given in a closed form, its evaluation contains a 
multi-dimensional integral over the A-body density matrix, 

p A (ri, r'i; r 2 , . . . , r A ) = ^ A ( r ii r 2, • • • , r A )^ A (r\, r 2 , . . . , r A ), (2.15) 

which is, of course, impossible for nuclear matter (A — > oo) . We previously carried out 
such a multi-dimensional integral for a light nucleus, 16 O, in a similar expression for 
the nuclear transparency [[2lJ, and examined various approximations, which could be 
used for heavier nuclei including nuclear matter. A possible approximation is to choose 
the A-body density matrix as 

A 

p A (r l , r\; r 2 , . . . , r A ) = p(r h r\) JJ p(r i )^(|r' 1 - Tj\) g(\r l - r^), (2.16) 

i=2 

where p(r, r') is the one-body density matrix with p(r) = p(r, r'), and the one-body 
density is normalized to one, i.e., J dr p(r) = 1. g(|rj — r 3 -|) describes both dynamical 
and statistical two-body correlation between the nucleon, 1, and nucleon, j. This 
approximation consists of including only the two-body correlation and of neglecting the 
so-called spectator effect |I2"|, p2|, which amounts to dropping two-body correlations 



among the nucleons, 2, . . . , A, and has been found accurate enough when the finite 
interaction-range of VNN( r ) is taken into account With the above A-body density 
matrix, the multi-dimensional integral over r 2 , . . . ,r A factorizes, i.e., 

A 

dr 2 ■ ■ ■ dr A exp{ — / dz" E ^NN(r"i - r,)} x p A (r 1 , r\; r 2 , . . . , r A ) 

V Jz 'i j=2 
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~ p(ri,r'i) 

x I J dr 2 p(r 2 ) gdr', - r 2 |) ^fln - r 2 |) exp{-^- jT dz'[ V^{v\ - r 2 )}j 
= P(ri, r'i) x[J dr 2 p(r 2 ) g{\v' x - r 2 |) g>(|ri - r 2 |) 

x (1 - 1 + exp{-- f 1 dz'l VM*"i ~ r 2 )})] A - 1 
~ p(r 1; r'i) x exp[-(A - 1) / dr 2 p(r 2 ) #(|r'i - r 2 |) g(\r 1 - r 2 |) 

x(l - exp{-- f 1 ^' V NN (r"i - r 2 )})], (2.17) 

where we have used the relation between the A-body and the one-body density matrices 
requiring 

J dr 2 p(r 2 ) g(\r\ - r 2 |) g(\n - r 2 |) ~ 1, (2.18) 

and assumed A to be large. 

Thus, we obtain the expression for the response function in the following form, 

S(u,q) = — f db 1 dz 1 dz[ p(ri,r'i) exp{i(p - q)(z 1 - z[)} (2.19) 
x [9( Zl - z[) exp{F(b i; Zl , z[)} + 6(z[ - z x ) exp{F*(b i; z[, Zl )}} , 
where F(hi \ Zi, z[) expresses the effect of the FSI (FSI-function), 

F(b 1 - Zl ,z' 1 ) EE -(A-l)JdT2 P (T 2 )g(\T' 1 -r 2 \)g(\Ti-T 2 \) 

ex P{"^ £ dz'{ V m (r'\ - r 2 )}^) , (2.20) 

and b'i = bi due to the eikonal approximation, eq. ([2.11|) . Note that we have not 
used the expansion in terms of the nucleon-nucleon potential, because the effect of 
the interaction with each nucleon cannot be assumed to be small. The expressions, 
eqs. ( [2.19 ) and (|2.20| ), are our main results of this work. 



One characteristic point of our formulation appears in this expression. Reflecting 
the fact that the Green's function for the struck nucleon, eq (|2.9|) , is off-diagonal with 
respect to r x and r'i, the correlation function appears as a product of g(\r'i — r 2 \)g(\ri — 
r 2 |). This feature is missing in the formulation based on the optical potential We 



will discuss the difference in the subsec. 3.2. 



2.3 Expressions for Nuclear Matter 

For later convenience, we here derive a simplified expression for the response function, 
eq. (|2.19|) , for nuclear matter. It becomes 

Y roc 

S(u,q) = — / dZ W(Z) exp{i(p-q)Z} (2.21) 

71 V J-oo 

x [6(Z) exp{F(Z)} + 9(-Z) exp{F*(-Z)}], 
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where Z = zi — z[, and V is the volume of the system. Since we normalize the one-body 
density to one, p(r) becomes po = 1/V. The FSI-function, eq. ( |2.20|) , is expressed as 

F(Z) = Pnm J dr 2 g(\r\ - r 2 |) flrflri - r 2 |) 

x fl - exp{-^- jf 1 dz'i V^"i ~ r 2 )}J , (2.22) 

where pnm = Ap is the nuclear matter density, and we use the value, pnm = A/V 
= 0.17 [fm~ 3 ]. Since both the FSI-function and the one-body density matrix become 
functions of the difference of the arguments for the uniform nuclear matter, we simply 
write F(Z) for the FSI-function and W(Z) for the one-body density matrix in eq. Q2.21|) , 

W(Z) = W(r-i f ) = p(T,T f ). (2.23) 

Using the fact that W(Z) is an even function, i.e., W{—Z) = W(Z), and introduc- 
ing a dimensionless function, w(Z), defined by 

w(Z) = VW{Z)=p(r,r')/p , (2.24) 

we further reduce the expression in the following way, 

y roc _ 

S(u,q) = — / dZW(Z) 

TTV JO 

x [exp{«Go - q)Z + F(Z)} + exp{-i(p - q)Z + F*{Z)}\ 

POO 

/ dZ w{Z) exp{Re F(Z)\ cos{(p - q)Z + Im F(Z)\, (2.25) 
Jo 

where we have decomposed F(Z) = Re F(Z) + ilm F(Z). This is the expression that 
we use for the numerical calculations for the inclusive cross sections in subsec. |4.2| . 

To compare our formulation with the formulation based on the optical potential 
in subsec. |3.2| . it is convenient to define the "convolution" function, ({Z). For this 
purpose we introduce the PWIA response function, S (u>,q), which is defined by 

1 f°° ~ 
So(u;,q) = — / dZ w{Z) exp{i(p — q)Z} 

TTV J-oo 

/oo „ 
dZ S {Z) exp{i{p-q)Z}, (2.26) 
-oo 

where S (Z) = w(Z)/(nv) is the Fourier transform of So(a>, q). Note that S (Z) is ac- 
tually a function of (p — q) with p related to to by eq. (P-12|) , due to our approximation of 
neglecting the residual nucleus excitation energies (the fixed-scatterer approximation). 
Using Sq(Z), we can write the response function as 

dZ S (Z) exp{i(p - q)Z] ((Z), (2.27) 

-oo 

where ({Z) is our "convolution" function defined by 

C(Z) = 9(Z) exp{F(Z)} + 9(-Z) exp{F*(-Z)}. (2.28) 

C(Z) includes all the information on the FSI. If there is no FSI, then ((Z) = 1, and 
S(u,q) coincides with S (u,q). 



2 

TTV JO 
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2.4 Zero-range Approximations 

In this subsection we discuss the zero-range approximations for the nucleon-nucleon 



potential, VnnM, for the FSI-function, eq. (|2.20|) , which are often employed in the 
literatures. Within our knowledge, the approximations have never been examined in 
the Glauber approach. Here we will study them in two stages within our framework, 
and reduce the numerical work greatly. 

The first one is to make the zero-range approximation only in the ^-direction (ZR1). 
This approximation amounts to replacing Vnn(i*) by Vj_(b) 6(z). The b-dependence of 
the last factor on the r.h.s. of eq. Q2.2Q) factorizes and we obtain 

l-exp{-~ r^^NN(r"i-r 2 )} ~ r(bi - b 2 ) 9{ Zl - z 2 ) 9(z 2 - z[), (2.29) 

V J zt 



where r 2 = (b 2 ,z 2 ). The FSI-function, eq. (|2.20| ), becomes 



~i 



F zm (bi, z 1 , z[) = -(A-l)jdb 2 J t dz 2 p{b 2 ,z 2 ) g{\r' 1 - r 2 |) g(\r 1 - r 2 |) 

xr(bi-ba), (2.30) 

where T(b) is the Fourier transform of the nucleon-nucleon scattering amplitude, and 
the definition is shown in eq. (|A.5|) in the Appendix. ZR1 should be a good approxi- 



mation for the finite-range interaction if \z\ — z[\ is large, because the step function in 
the ^-direction due to the zero-range approximation will be hard to discriminate from 
the smoothed one over the range of the interaction for large \z\ — z[\. 

Though the full calculation of the FSI-function for the finite-range interaction is 
probably difficult, the behavior of the small \z\ — z[\, where | dz'[ Vnn(i""i — r 2 )| is 
small, can be looked into by expanding the l.h.s. of eq. ( |2.29j ) in terms of the potential 
up to the first order. We write the expression ( [2.20|) for F FR (b!; zi, z[) Q as 



F PR (bi; zi, z[) = — - (A - 1) p dz'l [ dr 2 p(r 2 ) ^ NN (r"i - r 2 ) 

x^lr'i-r.D^dn-r.l). (2.31) 
This expression is useful to see how i*FR,(bi; Zi,z[) and that of ZR1 are different in the 



small \z\ — z[\ region. The numerical results are shown in subsec. 4J. 

The second one is to make the zero-range approximation in all the direction (ZR2). 
One cannot simply replace Vnn(i") by Vq S(r), because the appearance of the (^-function 
in the exponential is meaningless, which does not happen for ZR1. We have to be 
careful for determining the expression of T(b) in this approximation. From eq. ( |A.8|) , 
T(b) is written as 

r(b) = — £± exp{-zq-b}/(q). (2.32) 

l\p\ J {27T) Z 



1 FR stands for Finite Range. 
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The zero-range approximation in the transverse directions implies that /(q) is inde- 
pendent of q. Thus eq. ( |2.32| ) becomes 



r(b) 



2vr 

*|P| 
2tt 

*IpI 



/(0) <^(b). 



exp{ — zq ■ b} 



^2.33) 



Substituting this formula into eq. ( p.30[) , we obtain 

F Z R2(bi;^,4) = -(A-l)|^/(0) |dr 2i o(r 2 )^(|r , 1 -r 2 |)^(|r 1 -r 2 |) 



x5( 2 )(b!-b 2 ) e(z 1 -z 2 ) e(z 2 -z[) 

2n - . r z i 



(2.34) 

-(A - l)^, /(0) r dz 2 p(b 1; z 2 ) yflsj - z 2 \) g(\ Zl - z 2 \). 
The two approximations are numerically compared in the form of the FSI-function in 



subsec. EO and in the inclusive cross section in subsec. 4.2 
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3 Comparison with other Formulations 



In the following two subsections, we compare our formulation based on the Glauber 
approximation with the theory of Gersch, Rodriguez, and Smith (GRS theory) ||, || 
and a formulation with the optical potential [pTTJI . For the latter case, we show some 
numerical estimates in the next section. 



3.1 GRS Theory 

In this subsection, we point out the features of the GRS theory [§], J|, and compare 
them with those of our formulation by using R(uj, q) defined in eq. (|2.5|). 

R{u>, q) = (Ap(q) 1 6(q) \A), (3.1) 

UJ — ti A + IT] 

where O(q) (= exp{^q • ri}) is the hard scattering operator defined in eq.( |2.6|) , and 
Ha \A) = 0. The response function is obtained by taking its imaginary part as in 
eq. (|2.4|) . We decompose the full hamiltonian, Ha — K\ + Ha-i + Vi t A-i, in the same 
way as in eq. ( [2.7| ) of sec. |2]^. By expanding R(u>, q) in Vi^-i up to the first order, we 
compare the Glauber theory with the GRS theory. 

First, we consider the expansion of R(uj, q) in the GRS theory. Since the hard 
scattering operator, O(q), is the momentum-shift operator, it operates on the Green's 
function in eq. ( |3.1| ) to shift the momentum. The kinetic-energy operator of the struck 
nucleon, K\ = p 2 / (2m), is the only term affected in the hamiltonian, Ha, and is shifted 
as 

6t(q) K X 6(q) = K 1 + + (3.2) 

m 2m 

Substituting this expression into eq. (|3.1|) , we obtain the expansion in terms of Ha, 

*V,q) = (A\ ^ — VJ7^-^ 

uj — Ha — pi • q/m — q z / (2m) 

= (A\ : 7^ T77 AA) 

uj — p! • q/m — q 2 / (2m) 

+ (A\ ; ^ YT7^ ** A r YTT^\ A ) + " " " 

u> — pi • q/m — q z / (2m) uj — pi • q/m — q^/ (2m) 

= R (u,q) + R 1 (u,q) + ---. (3.3) 
Using Ha \A) = 0, we can rewrite R\(uj, q) as follows: 

- u; — pi ■ q/m — q z / (2m) 



1 


^ - pi 


q/m 


— q 2 / (2m) 




1 




u - pi 


q/m 


- q 2 /(2m) 




1 





( A \ = r 27?o v [#a, ; jtt^] l^> 

" •• uj — pi • q/m — q 2 /(2m) 

(-41 — : r — jt^tt [^m-i: — ; r — ^r^] l A >- 

— p! • q/m — q 2 / (2m) uj — pi • q/m — q^/ (2m) 
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Only V\ } a-\ does not commute with pi • q/m. 

Ro(u!,q) and Ri(u,q) in eq. ( |3.3|) gives F (q) and -Fi(q) of Rinat and Taragin 
T0|1 . They apply the GRS theory to the (e, e') inclusive reaction for the first time. In 



their work, after they obtain the expression of F (q) and -Fi(q), they exponentiate the 
potential term of -F\(q) to obtain their final expression. 

Next, we consider the expansion of R(u, q) in terms of Vi^-i m the Glauber theory. 
According to the usual perturbation theory we write 



u — Ha uj — K\ — Ha-i 

+ J—, V x , A -i + ■ • • . (3.5) 

to - Ki- H A -i to- K 1 - H A -i 

We apply the fixed-scatterer approximation and the eikonal approximation to the above 
expression. By the fixed-scatterer approximation we imply the replacement with Ha-i 
to Ea-i, where Ea-i is the average value of the excitation energy of the residual 
nucleus. By the eikonal approximation we imply 

(A\6\q) J——0(q)\A) = (A\ ^— 1 \A) 

u) — K\ — Ea-i uj — K i — px ■ q/m — q 2 /(2m) — Ea-\ 

- <A ; \ /{9 - \A). (3.6) 

to — p! • q/m — q z / (2m) — Ea-\ 

To obtain the second equality, we have put — 0, because the initial momentum of 

the struck nucleon is small compared to q. Corresponding to the expansion of eq. ( p.5[ ), 
R(u, q) of eq. (|3.1| ) is expanded as 

R(u,q) = R Q {u,q) + R 1 (u,q) + ---, (3.7) 

where Ro(u,q) and Ri(u,q) are written as 

Ro(u,q) = (A\ 1 - \A), (3.8) 

uj — pi • q/m — q z /(2m) — Ea-\ 

fli(w,q) = (A\ 1 V 1;A -i 

to — pi • q/m — q z j (2m) — Ea-i 

' \A). (3.9) 



a; — pi • q/m — q 2 / (2m) — Ea 

Now let us discuss the comparison of the GRS theory with the Glauber theory. For 
up to the first order expansion we can prove the following equation: 

R (u,q i ) + R l (u,q) = R (u, q) + R x {u, q) + A^ou, q) + A 2 Ro(u, q), (3.10) 

where AiRq(u, q) is the correction to the eikonal approximation defined by 

A 1 R (u;,q) = (A\- ^ U2 ^ilA), (3.11) 

{u> — pi • q/m — q £ l(2m)Y 
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and A2Rq(uj, q) is the correction to the fixed-scatterer approximation defined by 

1 



A 2 R (uj } q) 



{to — • q/m — q 2 /(2m)} 2 



(H A ^-E A ^)\A). (3.12) 



The proof of eq. ( |3.10| ) is as follows: The difference of the zeroth order in Vi )A -i is 
written as 

1 



R (v,q) - Ro(u,q) 



Adding the two correction terms, eqs. 
eq. (|3.13|) , we obtain 



E A ^\A). (3.13) 
{lu — pi • q/m — q 2 /(2m)\ 2 

|3. 1 1|) and (|3.12|), to the zeroth order expression, 



R {uj, q) - Rq(uj, q) + AxRoiuj, q) + A 2 R {uj, q) 

(^b ^ T~ 9//o MP (^i + ^a-i)|A) 



i?i(a;,q) 



p! ■ q/m — q 2 / (2m) } 2 
1 

— pi • q/m — q 2 /(2m)} 2 

- £i(w,q), 



Vl,A-l|^) 



(3.14) 



where we have used iY^ |A) =0. This completes the proof of eq. ( |3.10|) . 

The relation of the two formulations, eq. ( |3.10| ), implies that the expansion up to this 
order is exact in the GRS theory, while some corrections are needed to equate them in 
the Glauber theory. Thus one can realize that up to the first order in V\^ A -\ the GRS 
theory includes less approximation than the Glauber theory, though the corrections 
should be small in the high energy region. 

If we go beyond the first order, the situation will change. In the GRS theory it 
looks difficult to estimate higher order terms, while higher order terms can be summed 
up in the Glauber theory. This is the crucial difference, because the convergence of the 
series is not good from our discussion below eq. (|2.20|) in subsec. 22. Since the higher 
order terms are important, we should sum them up. Furthermore, the fixed-scatterer 
approximation and the eikonal approximation employed in the Glauber theory are 
established in high-energy regime. In this sense, we believe that the Glauber theory is 
superior to the GRS theory for descriptions of the high energy reactions. 



3.2 A Formulation with Optical Potential 

The optical-potential formulation has been used by various authors to treat the effects 
of the FSI in inclusive scattering (see e.g., Refs. 



17 , IS , 23 



'3)- Here we 

briefly review the treatment of the FSI and of the nuclear correlation in the formulation 
of Benhar, et al. ]ll[ . The main difference from our formulation will be seen to lie in 
the treatment of the nuclear correlation. 

In the optical-potential formulation the effect of FSI is included in the convolution 
form. The target nuclear tensor, W^ u (q, u), which corresponds to the response function, 
S(u>, q), in our approach, is expressed as 



/oo 
dJ F{u-uj')W^ IA {q,J-V{q)) 



(3.15) 
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where W^ ulA (q, uj) is the target nuclear tensor in the impulse approximation which 



gives the PWIA cross section corresponding to So(u;,q) in eq. ( |2.26|) . F(uj) is the 
convolution function expressed as 



F(uj - uj') 



2tt 



dt e i{ ^- J ^ e-wivtM 



-Re f°° dt e^-^* e- w ^* 

TT JO 

dt cos{(uj - uj')t} e- w(q ' t]t . 



TT JO 
j r 00 



TT JO 



(3.16) 



V(q) — iW(q,t) is an "optical potential". The t-dependence of the imaginary part is 
introduced to take account of the initial-state correlation between the struck nucleon 
and the other nucleons. This is thus not the usual optical potential to be used for 
the elastic scatterings and its foundation is not clear. The cu-dependence of V(q) and 
W(q, t) is assumed to be absent. If W(q, t) = 0, then F(u — uj') = 5{uj — uj'). 

The above expression, eq. ( |3.16| ), can be directly compared with our convolution 
function, eq. (|2.28 ), because both of them contain full information of FSI, and appear 
in the convolution form with the PWIA response function. The numerical comparison 



will be shown in subsec. fj~T . 

Let us write z = vt instead of t as a variable of W(q,z). Using the zero-range 
approximation for VnnOO, the authors of Ref. O] give 



W{q,z) 



^PNuv(q)a t N t l (q)- f ' dzf g(z'). 
2 z Jo 



(3.17) 



g(r) is the pair-distribution function, which corresponds to our correlation function, 
defined by [p5 



9( r ) = Pnn(pnm,^)/pnm- 



(3.18) 



Pnn(pnm, ?") is the average two-body density in nuclear matter at density Pnm- The 
average two-body density is defined as 



Pnn 



r 12 ) = T / PNN(ri,r 2 ) 



(3.19) 



with R 12 = (ri + r 2 )/2 and r 12 = ri — r 2 . For nuclear matter pNN( r i 2 ) is a function of 
|r 12 |, If there is no correlation, i.e., g(r) = 1, then 



W(q,z) = -p NM v(q) a™ at (q). 



(3.20) 



We cannot directly compare the correlation function defined in eq. ( |3.18j) with our 
correlation function defined in eq. (|2.16j) . However judging from how they appear in 
the term of the FSI effect, W(q, z) of eq. ( gig) and the FSI-function of eq. fl2~20D , it 
would be appropriate to identify their g(r) with our g 2 (r). 
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4 Numerical Results and Discussion 



In this section we will show our numerical results and the discussions. 



4.1 "Convolution" Function 

To calculate the convolution functions defined by eq. ( |2.28j ), we need a specific form of 
the correlation function, g(r). We take the correlation function of the form 

g(r) = 1 - Cl e-^, (4.1) 
where C\ = 0.84 [fm], and ai = 0.7 [fm] (Fig. 3), which roughly simulates that of 



Benhar, et al. (Fig. 6 of Ref. |l I]). In our figure we plot g 2 (r) instead of g(r) in order 
to compare it easily with their pair-distribution function, because our g 2 (r) corresponds 
to their g(r) as we discussed at the end of subsec. |372|. As one can see from those figures, 
the way how our correlation function approaches unity is slightly different from theirs. 
The difference in the curvature may cause a nontrivial effect on the cross section. In 
this case, we need another Gaussian term in eq. ( |4.1| ). Since what we would like to see 
is the difference coming from different manipulations for the FSI and the correlation, 
we believe that our choice of the correlation function causes no big problem. 

One more thing which we should comment on is the normalization of the wave 
function. For the case of a finite nucleus, we should be very careful for the normalization 



including the effect of the nuclear correlation Pq| . Fortunately the correction should 
be the order of 1/A, which can be neglected for nuclear matter. 

By using eq. ( |4.1|) , the approximate expressions for F(bi; z±, z[), eqs. ( |2.30|) and 
Q2.34Q , can be calculated explicitly for nuclear matter. The expression for ZR1, eq. ( |2.30|) , 
where the zero-range approximation is applied only in the z-direction, becomes 



F zm (Z) = FzRi(b i; (4.2) 



-Pnm J dh 2 J dz 2 g(\r' 1 -r 2 \) g(\r 1 -r 2 \)T(h 1 -h 2 

2«|p|7 l/af + 1/(47) V 2 a iJ 

+ 2/af + 1/(47) V=2 ^ I-23J { ^ U' rJ }] ' 

where Z — Z\ — z[ and (A — l)p ~ Pnm- Here T(z,p) is the incomplete Gamma 
function defined by 

r(z,p) = dt e _ V _1 , (Re z > 0). (4.3) 
Jp 

A useful formula related to the incomplete Gamma function is 

1 /I \ 

dxe~ x2 = - T (-,a 2 ) . (4.4) 
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The expression for ZR2 where the zero-range approximation is applied in the whole 
direction, eq. ( |2.34|) , becomes 

Fzr2(Z) = i ? Z R2(bi; z 1 , z[) 

f(0) Pnm / dz 2 g{\z x - z 2 \) gQz-t - z 2 \). 



*|P| 

/(0) pnm [Z - c x a x {\Ar- T ( ~, )} 



+c? ^ xp (-!i) { ^- r &5) }1 - (4 - 5 » 

Using the Gaussian series of the potential, eq. ( |A.12j ), in the Appendix, we write 
the expression for FR, eq. ( [2.31| ), as 

F FR (Z) = F FR {W,z h z[) (4.6) 
= -- Pnm / dz'l / dr 2 V NN (r'\ - r 2 ) g(\r\ - r 2 |) g(\rx - r 2 |) 

{(n/(4 7 ))(l/qf)}-V 2 /l (n/(4 7 ))(l/a 2 ) 2 V, 

Cl n/(4 7 ) + l/af W ^2'n/(47) + l/a? Z J 11 ' 

Note that this expression is valid only for small- Z (= Z\ — z[) region. From the 
discussion at the end of the Appendix, the expansion should be valid for \Z\ < 0.4 [fm] 
for the parameterization which we use below. 

To clarify the relation of ZR1 with FR we introduce one more expression, ZR1'. 
This has no physical implication, but is meaningful for obtaining some insight on the 
relation between ZR1 and FR. In order to obtain ZR1', we apply the expansion in 
terms of the potential to T(b) in eq. ( |4.2j ), 

T(b) = l-exp{ — / dz V NN {b,z')} 

V J— oo 

dz'V NN (b,z'). (4.7) 

i 

This approximation is similar to that of FR, but the integration range is different. 
Substituting this expression for eq. (|4.2|) , we obtain, 

F 2 m'{Z) = FzRi'Cbi;^,^) 

Jdh 2 J^ dz 2 g{\r\ - r 2 |) g(\r x - r 2 |) 

/oo 
dz'l VM*i ~ r 2 ) (4.8) 
-co 

Pnm Vtt 2^ V n \ ~ — ^Z - — ^- — — a x {^/n - T -, — } 

cjTT at ( Z 2 \ r n /l Z 2 \ . 

+ 2/a ? + ./(4 7 )71 eXP ™ { ^- T [r^ }] - 



V J-oo 



% i- r f zi 

— Pnm v 71 

v 



16 



The first term of this expression is the leading term for large-Z. This term is the same 
as that of FR of eq. (|4.6|) . Therefore, ZR1' and FR should show the same large- Z 
behavior. 

The numerical results for the convolution functions are shown in Figs. 4 a)-c). In 
Figs. 4a) and b) we plot 

Re C(Z) = exp{Re F(Z)} cos (Im F(Z)) , (4.9) 



because this form appears in eq. ( |2.25| ) for p = q, and it carries full information on 
FSI. In Fig. 4a) the results with the nuclear correlation are shown. The important 
finding here is that the slope of FR at Z = is the same as that of ZR1. This implies 
that the ZR1 can be a good approximation of the finite-range interaction case even for 
small-Z. Since we consider ZR1 as a good approximation of FR for large-Z, we expect 
it to be so for the whole range of Z. ZR2 shows a similar behavior to ZR1 in large- Z 
region, but in small- Z region ZR2 behaves differently from ZR1. Thus ZR2 cannot 
be expected to be generally applicable. ZR1' shows the same behavior as FR in the 
large- Z region, and slightly different in the small-Z. This is what we expected for the 
small- Z expansion and thus justifies our expectation that ZR1 is a good approximation 
for the whole region. 

In Fig. 4b) the results without the nuclear correlation are shown. ZR1 and ZR2 
degenerate as one sees from eqs. (|4.2|) and ( f4.5|) , and they show different behavior from 
FR even in small-Z. From those results, we realize that the nuclear correlation makes 
the zero-range approximations better. The nuclear correlation plays an important role 
in this context. 



In order to compare our formulation with that of the optical potential |TT 



we 



make the Fourier transform of our convolution function, ((Z), because their results 
were shown function of u, 

~ 1 r°° 

C(o;) = — / dte^((t) 

1 r°° 

= — dt (exp{+icut + F(t)} + exp{-icot + F*(-t)\) 
2n Jo 

1 r°° 

= - / dt exp{Re F(t)\ cos (wt + Im F(t)) . (4.10) 

TT Jo 



Here we use vt = Z. 

The numerical results are shown in Fig. 4c). In the figure we plot the Fourier 
transform of ZR1, ZR2, and ZR2 without the real part of Vnn(^)- This figure should 
be compared with Fig. 4 of Ref. [11]. The characteristic feature of our results of ZR1 
and ZR2 is the positive slope at u — 0. This behavior comes from the appearance of 
Im F(t) in the argument of cosine of eq. ( |4.10j) , and is the effect of the real part of 
Vnn(?")- To confirm this statement, we plot ZR2 without the real part of Vnn^)- The 
curve turns out to be monotonically decreasing. 

For reference, we plot ^op(^) defined by 

~ 1 r°° 

CopM = - dt cos(ut) e~ w ^ t)l . (4.11) 

7T JO 
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which corresponds to the case of Im F(t) = in eq. ( fOU| ), and is the same form 
as eq. ( |3.16|) , i.e., the correlated Glauber formulation of Ref. [JO]]. We use eq. (|3.17|) 



for W(q, t) with our parameter set, and g(r) in eq. (|3.17|) is replaced our correlation 
function, {g{r)} 2 of eq. ( [4.1| ). Therefore, the following expression of W(q,t) is not 



completely the same as that of Ref. [1 I 



W(q,t) = -Pnm^^T / dt'{g(vt)} 2 
Z t Jo 

= ip™,«'[i-c^ { ^-r(i,^)} 

+ c? * { ^-r(i,?!£W (4.i2) 

The g-dependence of W(q, t) is implicitly included in v and o"nn- As one can see from 
Fig. 4c), Cop( w ) shows a quite similar behavior to the result of ((uj) of ZR2 without 
Re Vnn except for the magnitude. 



4.2 Inclusive Cross Section 



We numerically calculate the cross sections of the (e, e') inclusive scattering off nuclear 
matter. Since we are interested in the treatment of the FSI and the nuclear correlation, 
other information, such as the one-body density matrix, is taken as an input. 

We determine the shape of the one-body density matrix, p(r, r'), from numerical 
results of the momentum distribution of a nucleon, iy(k). The one-body density matrix 



is related to the momentum distribution as [26] 



p(r, r') = Po J 



dk 



3 — ik-(r— r') 



W(k), 



(4.13) 



where po = 1/V. 

For the momentum distribution of the nuclear matter in eq. (|4.13|) , we use the 
numerical results of Refs. IE7I |2 



They use a realistic nuclear force to obtain the 



results. We just use their numbers by x-square fitting in terms of the following series: 



W(k) 



(2tt) 3 no(fc) 



(4.14) 



Ni 



(2tt) 3 \ e(k F -k){a + Y. 



n 9t j e jk2 / k F} + Q(k - k F ) ^ n e,j e" 



N 2 

E 



■jk/k 



where a = 7.156 x 1(T 2 , k F = 1.33 [fin -1 ], and k = 0.588 [fm -1 ]. iV\ = 15, and N 2 
= 10. Since the bare data are given by no(k) in Refs. |27| p^j , we show the relation 
of n (k) with our W(k) in eq. ( |4.14j ). We impose the condition that the Gaussian 
series and the exponential series are to be connected continuously at k = kp when we 
calculate the coefficients. The normalization is 



dk 



W(k) = 1. 



(4.15) 
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The resulting momentum distribution is shown in Fig. 5. For comparison, we plot the 
momentum distribution of the Fermi gas, W(k) = Wq 8(kF — k). The tail shows a 
reflection of the nuclear correlation. 

The numerical result for the one-body density matrix obtained from eq. ( |4.13| ) is 
shown in Fig. 6. For comparison, we plot that of the Fermi-gas model, which takes the 
form [EB1, 



P(r,r)=p ^-^, (4.16) 
where R = |r — r'|, and ji{x) = (sinx — xcosx)/x 2 is the spherical Bessel function. 



Putting the above ingredients into eq. ( |2.25 ), we obtain the electron inclusive cross 



sections. The numerical results are shown in Figs. 7 a)-c). For the electron- nucleon 
cross section we use the so-called Rosenbluth cross section |29|| , and average over the 
proton and the neutron, i.e., 

,da e N. _ 1 I da ep da en \ . . 

W el ' on = 2(^ + ^rJ- (417) 

The comparison of the calculated cross sections, ZR1 and ZR2, is shown in Fig. 7a). 
For comparison, we plot the results of PWIA and the case of no correlation. The 
PWIA result directly reflects the momentum distribution of a nucleon in the nuclear 
matter. Inclusion of only the FSI, i.e., g(r) = 1, broadens the cross section, but adding 
the nuclear correlation (ZR1 and ZR2) reduces the broadening. That is, the nuclear 
correlation reduces the effect of the FSI. Those are the general features of our results, 
and were also observed in Ref. [ IT| . By proceeding the zero-range approximation (ZR1 



— > ZR2), the cross section comes lower in the low-u; region. 

In Fig. 7b) we plot the cross sections for the cases of ZR1, ZR2, and no correlation 
with or without the real part of VnnO") to see the effect of the real part. As one can 
easily see from the figure, the effect is small for all cases. From our discussion at the 
end of the previous subsection, the real part affects the behavior of the FSI-function 
(see Fig. 4c)), but it does not affect the cross sections. We have thus confirmed the 
validity of neglecting the real part, which has often been presumed in the literature. 

In Fig. 7c) we plot our numerical result of the cross section of ZR1 (times 0.4) 



with the experimental data. From the discussion in subsec. |4.1| , the ZR1 includes the 
major part of the finite-range effect, but the cross section of ZR1 still overestimates the 
experimental cross section for low-u region, although this is a consistent treatment for 
the inclusive process. This would be due to our use of the on-shell eN cross section. 



As we mentioned in the beginning of subsec. |2.1| , we should use an off-shell eN cross 
section for the elementary process since the struck nucleon can become far off-the-mass- 
shell in the final state. Actually putting a binding effect by shifting the nucleon mass 
to the smaller one shifts the cross section to large-tu. Even this simple manipulation 
can make the agreement of the numerical results with the experimental data better, 
but this is only a prescription. Since one cannot determine the off-shell cross section 
without a dynamical model describing the internal structure of the nucleon and its 
interaction with the nuclear medium, we avoid to get involved in the off-shell problems 
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in this work. The problems are interesting and important, but we leave the study for 
our future work. 

In high-cj region, the experimental data highly exceeds our results, because the 
inelastic channel of the eN cross section takes part in the data. To describe the cross 
section in such a high-cj region is out of our scope of this work. 
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5 Summary and Conclusion 



We have given here a method to deal with the FSI in the high-energy (e, e') inclusive 
scattering. Since we are interested in the quasi-elastic region, the relevant degrees of 
freedom are the nucleonic ones. We have applied the Green's function method, and 
used the Glauber approximation for the FSI. An advantage of this formulation is that 
the finite-range effect of the nucleon-nucleon interaction and the nuclear correlation 
are included in a systematical way. 

Though the method allows us to express the FSI effect in a closed form, the full 
calculation still requires a large amount of numerical works. We have thus examined two 
kinds of zero-range approximations for the nucleon-nucleon interaction, which greatly 
simplify the integrals involved in the closed form. The one called ZR1, which uses 
the zero-range approximation only in the longitudinal direction, has been found to be 
accurate enough for the actual calculations. 

The method has then been compared with the treatments of the FSI in the other 
approaches to the response functions such as that of Gersch, Rodrigues and Smith, 



which has been applied to the (e, e') process in Ref. f]10| , and the optical-potential 
approach used in Ref. [|ll[ . 

Our principal objective is to propose a unified framework of calculating the inclusive 
(e, e') and the semi-inclusive (e, e'p) responses with the particular emphasis on the 
treatment of the FSI. The results of calculation in the present formalism (ZR1) for the 
(e, e') cross section off nuclear matter show strong effects of the FSI especially in the 



low energy transfer (u) region qualitatively similar to those observed in Ref. \TL\. The 
calculated cross section overestimates the experimental one by a factor of about 2 in 
the peak region and by larger factors in the low-to 1 region. This would be partly due to 
our simple choice of on-shell kinematics for the eN elastic cross section. We have not 
discussed the difficult problem of constructing the off-shell cross section, which would 
be required if we were to make a serious comparison with the observed (e, e') response. 
Instead, we have simply pointed out that, in the inclusive processes, the final nucleon 
is much more off-shell than the initial one. Further studies on the in-medium cross 
section including the off-shell kinematics are necessary to draw definitive conclusions 
on the FSI in connection with the color transparency. 

We should critically comment here the approximations which we used in this work. 
We have relied on several approximations, whose validation needs further and extensive 
work. One approximation is that the total disregard of E n of the residual nucleus, 



which we have used it to obtain the response function, eq. (|2.10|) . We believe it is 



good though we have no quantitative basis for it. The other approximation is to our 



choice of the form, eq. ( 2.16 ), for the A-body density matrix, which is claimed to 
represent a reasonable approximation for the calculation of the nuclear transparency. 
However it might be inadequate for the calculation of the low-energy-transfer region of 
the inclusive cross section. These approximations need to be investigated carefully for 
their quantitative validity. 
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APPENDIX 



A Nucleon-Nucleon Potential 

In this appendix, we explain how we construct the NN potential phenomenologically. 
This is an extension of the method introduced in Ref. M. 

We deal with a two-body scattering in free space, p is a momentum of the incident 
proton, p' is that of the outgoing proton, and the momentum transfer q = p — p'. 
Please do not confuse them with the notation of the main part of this paper. In the 
following we focus on a kinematic region where p • q ~ 0. 

The scattering amplitude is defined by 

/(q) = % M J d% T(b) exp{*( P - p') • b} (A.l) 
~ /(0) exp{- 7 q 2 }, (A.2) 
where b is the impact parameter defined by 



and 



b + z^-, (A.3) 
P 



/(0) = ^#|P|- (A.4) 



The last expression, eq. ( |A.2j ), is a phenomenological fit for the forward scattering 
amplitude. T(b) in eq. ( |A.1| ) is the scattering amplitude in the coordinate space, and 



it has the following eikonal expression 

F(b) = 1 - exp{i X (b)}. (A.5) 
x(b) is the profile function defined by 

1 r°° 

X(b) = — / dz'V^{h,z') (A.6) 

V J-oo 

= -log(l-r(b)). (A.7) 

where Vnn(f) is the nucleon-nucleon potential which we would like to obtain here. 

With eq. ( |A.1|) , we analytically obtain the expression of T(b) for the phenomeno- 
logical fit of /(q), eq. ( |A.2| ), by the Fourier transformation, 



i\p\ J (2tt) 2 

~ r(0)e" b2/47 , (A.8) 

where 

r( °» = 2^, m {A3) 
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Following the original paper of Glauber pj, we obtain the potential by using the 
Abel integral equation through eq. ( |A.6|) , 

vl d r^u x(b) 



^nnW = -~ [°°bdb 
7i r ar J r 



\Jb 2 — r 2 



vid r „„iog(i-r ( b)) 

7r r dr J r i y/b 2 — r 2 

-~ rdy-\og(l-T(y" + r')) 
ii r dr Jo % 



dy , VL-r^/^ . (A.10) 



277ri Jo 1 - T(0) e (f 2+r2 )/ 4 T' 



where y = a/6 2 — r 2 . We numerically calculate eq. ( |A.1C ), and obtain the potential. 
We use the following numbers at |p| = 2.0 [GeV/c] IpCfl : 

a^ al = 43.8 [mb], c = -0.14, 7 = 3.37 x 10" 6 [(MeV) -2 ], v = 0.905. (A.ll) 

The numerical results are shown in Fig. 8. 

For calculations of the response function, it is convenient to express the nucleon- 
nucleon potential in terms of the Gaussian series, 

No 

CW = E^" ,r!/47 . (A. 12) 

We determine the coefficients by \ 2 fit with the number of terms, Nq ~ 50. 

With the above parameterization, eq. ( |A.11| ), the absolute value of the following 
integral, 

-- p dz'iV^h^z'l) (A.13) 

= ~- E v n M ex P {-^ hj} - {v (—, ±A - r ( I, ±4) }, 

In P1 4 7 11 2 1 I, 2 47 / V 2 4 7 / 
stays less than unity for |zi — z[\ < 0.4 [fm] with bi = and z\ = 0. 
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Figure Captions 



Figure 1 

Kinematic notations in the (e, e') inclusive reaction. 
Figure 2 

Comparison of the off-shellness of the nucleon in the medium. 5m 2 /m 2 ^ is plotted 
as a function of E/m^. The solid line is the case (I), where V = 300 [MeV], S 
= -350 [MeV]. The dashed line is the case (II), where V = [MeV], S = -50 [MeV]. 
The dash-dotted line is the case (III), where V = —50 [MeV], S = [MeV]. The sum, 
S + V, is kept to be the same. 

Figure 3 

Nuclear correlation function, g 2 (r), of eq. ( f4.1|) as a function of r. 
Figure 4 a) 

Comparison of the two zero-range approximations for the nuclear matter with the nu- 
clear correlation. Real part of our "convolution" function, ((Z), of eq. (|2.28| ) is plotted 



as a function of Z (— Z\ — z[). The solid curve is ZR1, eq. (f4~2| ), and the dashed curve 
is ZR2, eq. (P~5|). For reference, we plot ZR1', eq. ([4.9|), by the dot-dashed curve, and 
FR, eq. ([Q>D, for a limited range of Z, by the dotted curve. 



Figure 4 b) 

Comparison of the two zero-range approximations for the nuclear matter without the 
nuclear correlation. The two are indistinguishable in the figure. The meaning of the 
curves is the same as Fig.4 a). 

Figure 4 c) 

Comparison of the two zero-range approximations for the nuclear matter with the nu- 
clear correlation. Fourier transform of our "convolution" function, ((Z), of eq. (|4 . X 0|) 



is plotted as a function of u. The solid curve is ZR1, eq. ( |4.2| ), and the dashed curve 
is ZR2, eq. ( |4.5|) . For reference, we plot ZR2 without the real part of VnnM by the 
dot-dashed curve, and the one based on the optical potential, (op(u) of eq. ( |4.11| ), by 
the dotted curve. 

Figure 5 

Momentum distributions of a nucleon in the nuclear matter as a function of k [fm~ ]. 
W{k) = (27r) 3 rio(k). The solid curve is the result of the fitting, eq. ([4.14]) , and the 
dotted line is the case of the Fermi gas. The crosses are the data points from Refs. [27 

u 



Figure 6 

One-body density matrix of the nuclear matter as a function of |r — r'|. The solid curve 
is the case including the nuclear correlation. The dashed curve is that of the Fermi 
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gas. 

Figure 7 a) 

The cross sections of the inclusive scattering off nuclear matter. The dotted curve is 
the case of PWIA. The dash-dotted curve is the case of only the FSI (no correlation 
effect). The solid curve and the dashed curve are the full calculation including both 
the FSI and the nuclear correlation with zero-range approximations. The solid curve is 
ZR1, the zero-range approximation in z-direction only. The dashed curve is ZR2, the 
zero-range approximation in the whole direction. 

Figure 7 b) 

Comparison of the cross sections for the cases of ZR1, ZR2, and no correlation with 
(the solid curve) or without the real part of VnnO") (the dashed curve). 

Figure 7 c) 

Comparison of the numerical result of the cross section and the experimental data |l] . 
The solid curve is the case of ZR1 x0.4. 

Figure 8 

Nucleon-nucleon potential, eq. (|A.1Q ), as a function of r for the case of |p| = 2.0 



[GeV/c], and 7 = 3.37 [(GeV/c) 2 ]. The dashed curve is the real part, and the solid 
curve is the imaginary part. 
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